> 



(N 



Analytic Approach to Perturbative QCD 

B.A. Magradze* 

A. Razmadze Mathematical Institute, Georgian Academy of Sciences, 

Tbilisi 380093 Georgia 



Abstract 

The two-loop invariant (running) coupling of QCD is written in terms of 
the Lambert W function. The analyticity structure of the coupling in the 
Q \ complex Q^-plane is established. The corresponding analytic coupling is 

' reconstructed via a dispersion relation. We also consider some other ap- 



proximations to the QCD /3-function, when the corresponding couplings 
are solved in terms of the Lambert function. The Landau gauge gluon 
propagator has been considered in the renormalization group invariant 
analytic approach (lAA). It is shown that there is a nonperturbative 
^ \ ambiguity in determination of the anomalous dimension function of the 

. gluon field. Several analytic solutions for the propagator at the one-loop 

order are constructed. Properties of the obtained analytical solutions 
are discussed. 



i!B ■ 1 Introduction 

{iLiI One of the most important tasks in quantum chromodynamics (QCD) is to find the 

^ ■ momentum dependence of the invariant (running) couphng as{Q'^)- For large Q^, the 

perturbative approximations of the couphng are rehable, since the theory is asymptotically 
. free. Outside of the asymptotic region, the perturbative invariant coupling is in fact large. 

\ Moreover, it has wrong analytical properties: the unphysical Landau singularities § 

\ appear at small space-like momenta. As a consequence, renormalization group (RG) 

improved expressions for perturbation theory (PT) approximations of physical quantities 
also have incorrect analytical properties. In particular, the electroweak current-current 
correlation functions parametrized by the running coupling in PT do not obey dispersion 
relations (DRs). 

DRs are fundamental for proving many important results in quantum field theories 
0) B il- There exists a vast theoretical literature for the application of DRs in deep- 
inelastic lepton-hadron scattering. This problem has been studied in rigorous manner on 
the basis of Jost-Lehmann-Dyson representation ||^. Other notable examples are the e+e~ 
annihilation ratio R^+e-, and the ratio of hadron to leptonic r-decay widths Rr- These 
time-like quantities can be related to the hadronic two-point correlation functions through 
DRs. These DRs provide a well-defined method for definition of the running coupling for 
time-like momentum. 

In the infrared region (IR), the investigation encounters great difficulties due to a fail- 
ure of PT. Several methods developed based on dispersion relation [|, 0, Borel summation 
and renormalons [0, ^, renormalization scheme (RS) choices and the Schwinger- Dyson 
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equations (SDEs) [|T0|. Common belief is that the unphysical singularities could be elim- 
inated when one adds higher order perturbative and nonperturbative corrections [H, 11 



The widespread supposition is that infrared properties of the exact effective coupling 
may lead to color confinement in QCD. Most of the suggestions of a confinement mecha- 
nism imply the assumption that the invariant coupling has a singularity at = 0. For 
instance, in the linear confinement picture, it is assumed that asiQ"^) ~ Q~^, ~^ 0. So 
the Landau pole moves out of the space-like region and causality is restored It was 
confirmed that this asymptotic behavior is consistent with the SDEs and Slavnov- Taylor 
identities of QCD 
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On the contrary, other authors have argued that the invariant coupling should be fi- 
nite in the IR region. It is assumed that the coupling "freezes" at low energies due to 
an infrared-stable fixed point of the theory. Such a possibility has been occurred in the 
K scheme [|15], |T^, where the corresponding /3-function in the three-loop approximation 
acquires an infrared fixed point. The similar conclusion has been achieved using Fade 



approximant methods or the Banks-Zaks expansion [0. Note that infrared "freez- 
ing" can not be in accord with the above mention linear confinement mechanism. However, 
there is an approach to confinement based on the BRST algebra and the RG methods (the 
so-called metric confinement) |jl9|, |20|. A linear quark- ant iquark potential is not implied 



within this framework. Instead, an approximately linear potential is supposed |21 



One possible resolution of the "Landau-ghost" problem has been suggested long time 
ago in Ref. [|^]. In Ref.[^, the idea to combine the RG invariance and the Q^-analyticity 
(in the context of quantum electrodynamics (QED)) have been proposed. A RG invariant 
analytic version of FT has been elaborated in this work. Recently, this method has 
been successfully applied in QCD. It was shown, that the analytic running coupling is 
stable for the whole interval of momentum, and has an universal infrared limit at = 
The "analyticized" perturbative approximations, as opposed to standard FT, exhibit 
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reduced RS dependence and the results are not sensitive to higher- loop corrections [23 



These features of RG invariant analytic approach (lAA) have been demonstrated in recent 
calculations of Re+e- |24| and [^. The calculations have been performed to next- 
to-next-to-leading order. Other improvement of the analytic approach is that it provides 
a self-consistent determination of the running coupling in the time-like region [Q. In 
2^ the Bjorken and the Gross-Llewellyn Smith sum rules in lAA are considered. 

For most recent 



Refs. 



The Landau gauge gluon propagator in lAA has been studied in Ref. 1 33 



applications of the analytic approach to QCD see the recent review |31 



In Sec. 2 we consider some technical aspects of the RG lAA. The two-loop running 
coupling is solved explicitly in terms of the Lambert W function. This allows us to in- 
vestigate analytic structure of the coupling in the complex plane. The corresponding 
analytic solution is reconstructed via the dispersion relation. To discuss the RS depen- 
dence, we compare the exact expressions for the causal coupling in the two different 
renormalization schemes. In Sec. 3 the structure of the /3-function in AFT is studied and 
other approaches to the Landau pole problem are discussed. The special (nonperturba- 
tive) model /^-function is investigated. In this case, we also obtain the explicit solution for 
the coupling in terms of the Lambert W function and show that the solution is consistent 



2 



with causality. In Sec. 4 we apply the analytic approach to the gluon propagator in the 
Landau gauge. Several ways of restoring analyticity for the propagator are considered. 
The ambiguity in causal "analyticization procedure" is discussed Sec. 5 contains the 
concluding remarks. 



2 Lambert's W function in lAA 

The running coupling of massless QCD satisfies the differential equation 

Q'^a.{Q')=P{as{Q')), (1) 

with the initial condition 

= as, (2) 

2 

here = is the renormalized coupling constant, and /i is the renormalization 

point. We assume that the non-trivial (3 function exists with the perturbative asymptotic 
expansion in powers of as 

P{as) ^ -^{asY - T^ias)' (3) 



Att' ' (47r; 



the first two coefficients of the formal power series are independent of the chosen renor- 
malization conditions. Their values are 

= (ll - ^iV/) , A = (l02 - ^iV^) , (4) 

with Nf being the number of quark flavors. The solution of Eq. (|I]), which satisfies the 
initial condition has the form 



VAV (3Qas{Q^) /3o V (3ias{Q^)) 



where 



and denotes the two-loop /5-function. The QCD scale parameter A here is different 

from the conventional one A;^, so that A = (6)^~^Ajg^ with b = ^ |36|. To obtain 

explicit expressions for asiQ"^) as a function of Q^, beyond the one-loop order, one must 
solve the transcendental equation (||). 

"We use the notation = —q^, > corresponds to a spacelike momentuni transfer. 
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Let a^'^\Q'^) be nth-order perturbative approximation to as(Q^) determined implicitly 
by Eq. (^. Then the corresponding analytic running coupling can be defined via the 
Kallen-Lehmann integral ^ 



the spectral function p(")(a) is the discontinuity of the "initial" expression q;'^'^-*(Q^) along 
the negative axis: p'^'^\a) = Ima^'^\—a — zO). The remarkable result is that the 



limiting value aa„(0) is universal and determined only by one-loop calculation |22 

«i"J(0) = a«(0) = ^. 

Po 

By using Cauchy's theorem Eq. (^) can be rewritten as follows 

ai:)(Q^) = a(")(Q2) + 0(")(Q2). (7) 

Here the "nonperturbative" term 9'^'^'^ (which compensates the unphysical contributions 
of a^"-*) comes from the ghost cut 

TT Jkl [a + Q'^ — zO) 

where < 0, to first and second orders kj^ = — A^. 

Let us now consider the two-loop running coupling in more detail. It is more convenient 
to introduce the quantity 

a(x) = ga(Q2), where x = ^. (9) 



To second order Eq. (R) reads 



^ Mn(l + ^— J— I =lna:, (10) 



where we have denoted ^ = ^ = |f for = 3). In what follows we shall consider 
the phenomenologically interesting case when Nj < 8.05 (0 < 6 < 1). Then the two-loop 
running coupling has the ghost singularity on the positive x-axis. The transcendental 
equation (|l^) can be solved by the iteration method. The solution given by one iteration 

aSHx) = , : ^. (11) 



Inx + 61n (l + iln 



X 



''Historically, for the first time the 2-loop iterative solution has been proposed in Ref. |Q in the 
context of QED. 
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Consider function ( [Tl| ) in the whole complex x-plane. It has the Landau pole at x = 1, 
the logarithmic branch point at x = exp{—b) and standard branch point at x = 0. The 
corresponding analytic coupling is given by the Kallen-Lehmann integral (see Refs.[p^2|) 

= - r = i r -^^m. (12) 



TT Jo S + X TT J~oo e* + X 

where p^\t) ^ p^\e') ^ p'i\e'A') = with 



Li = TT + 6arccos(^-^), L2 = t + bln(^), r = J (b + tV + tt^. 
r * 

One can investigate Eq.(|TU|) using the lagrange inversion formula. Then one finds that 
the solution has second order branch point at x = 1: a'^'^\x) ~ , ^ . Hence, approxi- 

a/26(x-1) 

mation ( pT]) violates the analytical properties of a'^'^\x) near the point x=l. In fact, the 
transcendental Eq. (0) is exactly solvable ^ 0. The solution is 



a(2)(x) = -i— i— : u{x) = W{(), (13) 
1 + uj[x) 

where 

1 1 / ln(x) \ , , 

C = --x"^ = exp ( ^-l + mj, (14) 



and W{Q denotes the Lambert W function This is the multivalued inverse of 



C = PF(C)expPF(C). 

The branches of W are denoted ^^(C), = 0, ±1, .... A detailed review of properties 
and applications of this special function can be found in |^S[. We shall define branches 
of W following this work. Then the branch cuts are chosen as {C : —00 < C ^ ~~} ^"^^ 
{C : -00 < C < 0} g The branch W_i{Q satisfies iy_i(C) < -1 for ~e~^ < C < 0. 
Now our task is to choose the suitable branches of W in Eq. (p^). The solution a*^2-*(x) is 
supposed to have the following properties: 

i. It would be an analytic function in the cut x-plane with cuts {x : —00 < x < 0} and 
{x : —00 < X < 1}. 

ii. It is real and positive for a real positive x > 1. 

iii. a*^2-*(x) ~l/lnxifx^oo along all directions in the complex x-plane. 

It is obvious, that for real positive x > 1 the relevant branch is uj{x) = W^i{Q. We 
look for the analytical continuation of uj{x) to the whole complex x plane. To perform 
the analytical continuation, one has to take into account the ranges of the branches of 
w = W{C) 0. Let us consider the cases < Nf < 6.19 (0.5 < 6 < 1) and 6.19 <Nf<8 
{0 < b < 0.5) separately: 

'^These notations and definitions of branches of W are different from that used in Refs.lR4| 
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1. < Nf < 6.19. Evidently, on the upper half-plane {Im{x) > 0) the solution is 
determined by W-i{(). On the lower half-plane the branch should be changed. Indeed, 
we have crossed the cut in the complex (-plane. Accordingly, the replacement arg((^) =^ 
arg(C) — 2tt is implied. We conclude that here the relevant branch is Wi{(). With this 
choice, the solution is continuous on the line 1 < x < oo. Thus we obtain 

, ,(^\ = I ^-idCle'^'O : ^1 = - i arg(x) if < arg(x) < vr. 

1 Wi{\(\e'^-') : = -vr - i arg(x) if -vr < arg(x) < 0. ^ ^ 

Here |C| = ^kl"'^. 

2. 6.19 < Nf < 8. In this case, we separate the upper half-plane, (arg(x) G (0,7r)), 
into sectors 

2{n - l)bn < arg(x) < mm{2nbiT, vr), 1 <n< rimax, ^ < n^ax < ^ + 1, 

(for Nf = 8, Umax = 25). An image of a sector is a sheet of the Riemannian surface of (, 
which maps onto the range of a branch of uj{x) = W{(). The analytical continuation to 
the n-th sector can be written as 

uj{x) = W_n{\C\e''^-"), where v?-„ = (2ra - l)n - I arg{x). (16) 

To obtain the solution on the lower half-plane, (arg(x) G (— vr, 0)), one must use the 

formula W„(C) = W^;:(Oil. 

Now we can reconstruct the corresponding analytic coupling a^^^i^)- the case 



Table 1: The dependence of the analytic running couplings q;^^^(Q^), ait^niQ"^) 
<-(Q2) for Nj 
= l.mGeV 



om'^iQ ) for Nf = 3. Here, we have used the reference value of as^M^ ) = 0.36 for 



Q2 GeV^ 


ai'„^(Q2) 




«ar(Q') 




Q2 GeV^ 


airi(g2) 









1.396 


1.396 


1.396 




5 


0.330 


0.330 


0.330 


0.001 


0.985 


0.969 


0.980 




6 


0.319 


0.319 


0.319 


0.01 


0.846 


0.831 


0.838 




7 


0.310 


0.310 


0.310 


0.05 


0.716 


0.706 


0.710 




8 


0.302 


0.302 


0.303 


0.1 


0.654 


0.647 


0.649 




9 


0.296 


0.296 


0.296 


0.5 


0.507 


0.505 


0.505 




10 


0.290 


0.290 


0.291 


1 


0.447 


0.446 


0.446 




20 


0.256 


0.257 


0.257 


2 


0.393 


0.392 


0.392 




30 


0.239 


0.239 


0.240 


3 


0.364 


0.364 


0.364 




40 


0.228 


0.228 


0.229 


4 


0.344 


0.344 


0.344 




50 


0.220 


0.220 


0.221 



< A^/ < 6, with Eqs. (D, (|T3D and (|l|), we obtain 



1 



n Jo s + X vr 



(e* + x) 



(17) 
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Table 2: Numerical results for A and A{2GeV) with exact and iterative solutions. 





0.34 


0.36 


0.38 




0.606 


0.706 


0.815 


A(2)(GeV) 


0.665 


0.772 


0.889 


Aspec(GeV) 


0.644 


0.750 


0.865 


A\f\2GeV) 


0.476 


0.499 


0.522 


A(-^\2GeV) 


0.479 


0.502 


0.525 


Aspeci2GeV) 


0.477 


0.500 


0.523 



where p'-^^t) = p^^\e^) = p^'^\e^K^) and 

p(2)(t) = -\im {j^wkm) ' = (-^ - 1 + ^(i - 1)^) • 

Equivalently, using the Cauchy formula we can rewrite Eq. ([T7|), in the form 

1 /•! 1 

TT Jo S — X 

In particular, for < x < 1, we have 

1 /-i 1 



a(^2j (3;) = ^(2) {x)-^£ 73-;^^^ (s + zO) 



Re{a'^^\x + iO)}-PV- f Im{a^^\s + iO)}ds, 

TC Jo S — X 

where PV denotes the principal value of the integral, and 

a^^\x + ^0) = -\ ^- - 



Numerical results for the exact two-loop coupling (|T^) as well as for the iterative solution 
(|12D are summarized in the Table 1. The relative error for the iterative solution ([T2|) , is 
less then 1.8 % for the considered interval. This gives 8% error in the value of the QCD 
scale parameter A^^^ for Nf = 3 (see Table 2). Following Refs.[^] we may use the average 



1 fQ _ 

MQ) = / ^s{p )dp 
Q Jo 



Q 

for comparison. For Q = 2 GeV, we find that iterative solution (|12|) gives answer for 
A{Q) to much better than 0.5% accuracy. The reference values of aan{Mr'^) are taken as: 
aa„(Mr^) = 0.36 ± 0.02 with Mr = 1.777 GeV. The instructive example is given by the 



/3-function of the special RS 
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in this scheme the inverse /3-function contains only two terms (for application of this RS 
see.[3£]). Note that expression (|l^) is beyond the formal framework of PT where only 
finite order polynomials in the coupling are allowed 0. The RG equation (|I|) with (^) 
yields the implicit solution for the coupling 



(6a(^))^exp(^): x = f , b = ^, (20) 

here A we define following Ref.||36|. An inversion of (|20|), can be written in terms of the 
Lambert W function 

('('^) = -b^)- Mx)=W{z), z = ~x-i. (21) 

For X > e^, {b > 0, for Nf < 8) the physical branch is W-i{z). This branch is real and 
yields the correct ultraviolet behavior for the a{x). We see that a{x) has the branch point 
(the unphysical singularity) at x = e^, and is finite at this point 



X 



^ \ 1 1 2(x-e*') b 

Let us now perform the analytic continuation in the complex x-plane. For < Nf < 6 
we obtain 

(r.] = S W^-i(k|e'^+) : ^+ = TT - i arg(a;) if < arg(s) < vr. 

\ Wi{\z\e'^-) : ^_ = -n- iarg(a;) if -vr < arg(x) < 0. ^ ^ 

The corresponding causal coupling is defined via the DR 

1 f°° Pl{s) 1 6* 

aan{x) = - — ■ — ds = - — rPi{t)dt, (23) 

71 Jo S + X TT J-oo [e^ + X) 



where 



pi{t) = pi{e*) = -\lm ^ 



Let us compare the exact results ([T7|) and (^3]) (note that the exact two loop expression 
(|T^ is the coupling of the 't Hooft scheme). Numerical results are summarized in the 
Table 1. The difference between the solutions in the infrared region is less then 1%. 

There are several reasons to believe that the application of the Lambert function may 
be useful from practical point of view: 

1) Of course, the transcendental equation (^ can be solved in the complex domain 
using numerical methods ||2^, ^ . However, this implies a preliminary analytical investi- 
gation: It is necessary to determine the analyticity structure of the implicit function (the 
relevant branch should be chosen). The application of the Lambert W function facilitates 
this task. 

2) The explicit inverted 2- and 3-loop solutions of RG equation may also be useful in 



studies of renormalon properties of gauge theories |51 



Note added. While preparing this manuscript I became aware of the paper which 
has some overlap with the present work. The authors of |^ have obtained explicit 
expressions for the coupling (at the 2-loop and 3-loop orders) in terms of the Lambert W 
function. 
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3 The /3-function in the lAA 

It is instructive to consider the /3-function. From Eqs. (||), (|]) and (j^) we have 

(3t\c.s) = -- r r-^,P^''\^)da = -- r -T-%-,ds, (24) 
TT JO [a + fi^)^ IT Jo y[s+-)'^ 



y 

where p^'^\s) = p^^\a), s = and y denotes the dimensionless "nonperturbative" 
variable 

y = exp(-0) = — , (25) 

for as — > 0, y ~ exp(— -g^) . The spectral representation (p^) defines P{as) = (3{y) as an 
analytic function in the cut-y plane. The cut is along the negative j/-axis. Another useful 
representation for the /3-function follows from Eqs. (|^) and (|^) 

= ^ = Pi\o^^''\^^')) + = &{<^s - ^("H/^^)) + ^"^"H/i^), (26) 

where denotes the nth-order perturbative approximation for the /^-function 



/3ir(«) = -E 



.=0 (4^) 



{~^) I — \ " ^ k+2 



k+1 



a 



and 9^'^\pP') = p? — ■ -^l- (H) "^^^ be rewritten as follows 

Pan{0ls) = Ppticts) + Pupils, y). (27) 

Here, the nonperturbative piece, denotes the generalized power expansion for as 

P^J{c^s,y) = jZBt\y){as)\ (28) 

fe=0 

where 

Bo = e^ 51 = -^, B, = i^'-^ + j^ for fc>2. (29) 

Here, for convenience, we use condensed notations suppressing the superscript (n) and 
arguments to the expressions above {0 = O'^'^^jj?) ets.). From the spectral representation 
(H) and Eq. (^) we see that the coefficients Bkiy) are analytical functions of y in the 
cut |/-plane. The cut is along the real positive y-axis y > > . In particular, the 
coefficients are regular functions in the neighborhood of ?/ = 0. For as — > 0, the coefficients 
vanish exponentially Bk{y) y exp(— -g^). 

A detailed study of the one-loop /^-function can be found in Ref. It was shown 
that the one- loop /^-function obeys the symmetry property P{as) = ~ c^s), and it 



has the second order zero at = To second order the symmetry property for the 
/3-function is not occurred. 

Notice that in the RG lAA the nonperturbative exp{—-^)-tjpe terms play an essential 
role for restoring of analyticity. The reason is that one starts with a finite order pertur- 
bative approximation to the /3-function. Another point of view has been advocated in 
Refs. 0. It was shown that a perturbative series for the /3-function, in the specific class 
of schemes, can be summed. So that the resulting running coupling has causal structure. 
The nonperturbative terms are not added, but the freedom of choosing a RS for an infi- 
nite series was used instead [§. The connection between that approach and lAA is still 
unclear. Indeed, they offer two different causal approximations for the running coupling. 



To illustrate we consider the model /3-function introduced in Ref. fi^] 



m = -^n-- h = fi, m=^^Pi<ys), (30) 

which yields a linearly rising potential and asymptotic freedom. The RG equation with 
(|30D can be solved, the solution is 

MS) = ^^^^;^: A^ = Vexp(-^), (31) 

here Wo(z) is principal branch of W{z) Quite remarkably, (|3TD satisfies the cut-plane 
analyticity. Indeed, Wo{z) is analytic at z = {Wo{z) = z — z"^...) and its branch cut is 
{z : —00 < z < —1 / e} . A generalization of this model has been studied in Ref. ||14|| . 

Generally, the running coupling is a gauge- and scheme- dependent quantity (it is non- 
physical object in this sense 0). Presumably, more useful quantities, in the infrared 
region, are the scheme-independent effective charge (the combination of a vertex and 
propagators) |^) or the running couplings of physical schemes |T^, Most satis- 

factory approach can be based on the notion of scheme- and gauge-independent universal 
charge (the analogous of the QED effective charge) Q]. Recently, in Ref. [^] the universal 
one-loop QCD coupling has been constructed using the pinch technique. 

It is clear that lAA does not exhaust all nonperturbative effects of the invariant cou- 
pling because its origin is PT. In Refs. ^ Q more general frameworks are considered. 
It is assumed that a nonperturbative universal QCD coupling exists which satisfies the 
DR. 

We can write the nonperturbative universal coupling as a decomposition 
The infrared and ultraviolet pieces of the coupling satisfy the DRs 

-ifAQ') = ^ "^f ''^^-^ -7fm = J , ^od-- (32) 



The momentum K{as, n) divides the infrared and ultraviolet regions. It may be chosen in 
a RG-invariant fashion. In Refs. it was shown that this nonperturbative scale actually 
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exists in QCD if Nf < 9. In the weak coupling limit the scale vanishes exponentially, 

47r 

K(as,fJ,) ~ A ~ e ''o"» as as — > 0. The APT approach corresponds to the simplest 
possibility: it is assumed that K = 0, hence a'J^f j^Q"^) is ignored. In addition, the weight 
function p^^{ij) is approximated by PT. 

One can, in principle, examine the consistency of lAA with a truncated set of SDEs. 
The authors of Refs. |Q have given reasons to believe that alfj-^Q"^) cannot be zero, and 
an infrared regular QCD effective charge cannot be consistent with the SDE for the gluon 
propagator. Ghost-free axial gauge has been chosen in this work. On the contrary, the 
authors of Ref . PS| have found that the infrared- regular nonperturbative effective charge 
is consistent with the truncated SDEs in the Landau gauge. However, it was shown that 
the corresponding gluon propagator has not causal structure. We see that situation is not 
completely clear. 



4 The gluon propagator in lAA 



In general the gluon propagator is not observable. But, this fact does not imply that it 
does not contain physics. One can relate the propagator to gauge invariant quantities, for 
example the Wilson loop [^, the gluon condensate or the gluon distribution function 



of hadrons Although the propagator is known to be a gauge variant quantity it 
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In the standard 



contains an important information about the infrared region 

It has been proved that the propagators in QCD obey the DRs |^ 
RG improved PT the DRs are violated Therefore it is desirable to develop the 

analytical approach for calculating of the propagators. It was pointed out in Ref. p8 
that there is an ambiguity due to "noncommutativity" of "analyticization" with some 
elements of the RG algorithm. For this reason, in the case of Green's functions and 
observables, several different versions of "analyticization" are possible. To illustrate we 
shall consider the gluon propagator in the Landau gauge 



1. 



(33) 



We have assumed that quarks are massless and normalized the propagator at the Eu- 
clidean point = fi'^. In the Landau gauge, invariance under the RG leads 



d 



2 ' 



exp 



dx 



(34) 



The formal power series for the anomalous dimension function 7^, (in this gauge) is given 
by 

lv{as) = -(^«s + ^as' + ...), 70 = i(13 - |iV^). (35) 
A detailed investigation of the gluon propagator in the framework of analyticity and 



asymptotic freedom has been undertaken in Refs. [ 48 1. It has been shown that the propa- 
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gator amplitude satisfies an unsubtracted DR 



(36) 



For a limited number of flavors Nf < 9, in the Landau gauge, the weight function 
obeys the super convergence relation 



Pv{a, as,p)da = 0. 



(37) 



This relation is considered as a sufficient condition for color confinement in the approach 
to a confinement based on the BRST algebra and the RG [jl9|, |2^. The consequence of 
( |37D is that there is a renormalization invariant point K'^{as, p) such that is a negative 
measure for a > 

Let us consider few different scenarios of restoring analyticity for the gluon propagator: 
(a) The Analytic Perturbation Theory (APT) possibility. This specific recipe for 
"analyticization" of an observable has been introduced and elaborated in Refs. |2^, Q . 
Here, instead of the power perturbation series, an amplitude is presented in a form of an 
asymptotic expansion of a more general form, the expansion over an asymptotic set of 
functions [a"(a;)]a„, the "n-th power of a(x) analyticized as a whole" [^. In the APT 
approach, the drastic reduction of loop and RS sensitivity for several observables has been 
found. 

In this case the starting point is formula (0). In the one- loop approximation (0) 
yields 

p 



DW(^ a) = — 



47r 



70 



{31 



Since the ratio ^ is fractional number, expression (^) has the ghost cut < < A^ 
The discontinuity p^^ of D^^^ along the negative axis can be written as |^ 



n'^^U^) 1 70 

pi'Ks) = -cMs)-{R{s))^osm 



To 
(3o 



$(5) 



(39) 



70 

where s = = (^) , R{s) = (In^ s + vr^)"! and 



$(.) 



arcsin(7ri?(s)) if s > 1 

vr — arcsin(7ri?(s)) if < s < 1. 



To construct the corresponding analytic expression one has to substitute (|39D in formula 
(|36D. Note that the weight function (^) has a nonintegrable singularity at a = 0. So 
that the spectral representation ( p6D with the weight function ( p9| ) diverges. Nevertheless, 
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in sense of the distribution theory this problem may be solved Using the method of 
Refs. we find the following subtracted DR for the "analyticized" amplitude 



where 



ci = ^ J p^J,legis)ds; p^^j^g = pi^) - p^^li^g] p^^lingis) = —^T^^ sin(p7r), (41) 



with p = ^ = 2{33-2Nf) (0 < p < 1 if A^/ < 9), and pi^\s) is given by (|39D . Instead of the 
super convergence relation (|37| ) now we have the following relation 



jp^'lg{s)ds + j-p^\s)ds = Q. (42) 

1 

However, there is another way to handle the above considered infrared problem. The 
point is, that for Nj < 9 one can write an unsubtracted DR for the dimensionless structure 
d, d = Q^Di^Q"^) This allows us to construct the analytic expression for d. The 

corresponding "analyticized" amplitude Dan is then determined by 

n(i) /n2^ _ ^inC^') _ ^ 7 P^^hs)ds 

^''-'^^ ^ ~ ~qhI TTq^ ^^^^ 

where Q'^ = jj, Pvii^) = ~^P^vK^) > with p^^-'(s) given by Eg. (|39D . One can verify 
that expression (^) has the correct "abelian limit" . Indeed, we have 



70-/30 ' vrJ s + g2 Vln"s + 7rV ' " ' \\nQ^ 1-Q 



Using the relation (^2]) we verify that the analytic expression (^3]) follows from the repre- 
sentation (|iD|). 

Note that ( ^31) does not possess a pole at = and it has a more strong singularity 
then the corresponding amplitude for the free field. By direct calculation we find 

More serious difficulty is that there is not proper connection between the analytic prop- 
agator (^) and the analytic running coupling as it follows from basic RG relations . 
Indeed, the product of a vertex and appropriate powers of propagators should form an 
invariant charge. 

''The infrared behavior of spectral representations has been studied in Refs. |4q|. 
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(b) One can substitute the exphcit expressions aaniQ"^) in formula p^). The simple 
possibility is to use the perturbative expression (|35| ) for 7„. This choice has been proposed 
in Ref. Then in the one- loop approximation formula ( |5^ gives the causal expression 



g^ing^j 



70 
/3o 



(44) 



where Q = j^, and 0^2 is the normalization constant. As in the case (a), we see that (0) 
does not satisfy the unsubtracted DR (^). The lack of this solution is that it does not 
reproduce the "abelian limit" for 70 — >■ /3o- 

(c) We may introduce RG improved power asymptotic expansion differing from the 
usual one by substitution aperi.(^) =^ o-an.i^) only. With this choice, in the one-loop 
order, we get 



where c„ 



7q _ 

^^^^0 andQ2 = ^. It 

is easy to convince that (^) has correct analytical 

properties. Indeed, a2^{Q'^) does not vanish in the finite part of the complex Q^-plane. 
Moreover, (^51) satisfies the unsubtracted DR (see Refs. [^, 0). For the RG invariant 
scale K{g,fi^), the solution (|45|) yields trivial value K = 0. In a similar way, one can also 
obtain causal approximations to D{Q'^) at higher orders, 
(d) Consider the one-loop gluon self energy 



a 



(1) 



20 



2a 



a. 



Q^'ninQ2 + i_g2 



(45) 



^«(Q^a,) = ^a,lnQ^ 
Po 



(46) 



here 70 is given by (|35|) , 



£«s and 



We can introduce the Dyson series in 



n^^) for the propagator. Summing this series we get 

1 



l + n(i)(Q2,a, 



-K In y 



(47) 



here k = ^ and instead of we have introduced the variable y, y = e~°^. The approx- 
imation (^7]) does not satisfy RG invariance and has the ghost pole. We may define the 
corresponding "analyticized" amplitude 0] 



d 



Ds.an 



{Q\y) 



1 



CDs{y) 



1 



y 



ln(Q2^-) g2 



(45 



in order to preserve the normalization (see (|33|) ) we have introduced in (|^) the necessary 
factor CDsiy)- The effect of this procedure is that the gluon self-energy receives a pure 
nonperturbative contribution. The modified gluon self-energy function is then given by 



CL(0^z/) = ^«(0^-(ln^/)-^) + ^, 



(1) 

nonpert 



{Q\y), 



(49) 



14 



here 

with n(i)(Q2) = n(i)(g2, -(lny)~^) and G{Q^) = y''{Q^ - y'')-\ We see that the expres- 
sion (^91) is free from the ghost singularity. For Nf < 9 (k > 0), in the weak-couphng 
hmit, the nonperturbative part of the self energy vanishes exponentially: ^nonpert.iQ'^ ^ V) ~ 
0(e~"^) 0. The modified one loop amplitude is then given by 

rfw(Q^^/) = l-^;;L(Q^^/). (5i) 

Now, the perturbative definition of y, y = e~'^ , is not valid. Instead, we shall accept 
the modified relation y = ^ = e~'^^°"'^ (see Refs.[|], ^). Thus, to first order, (f){as) is 
determined by 



1 

+ 



We can calculate the modified anomalous dimension for the gluon field 

d d 
7mL(as) = 7mod(0) = Jim Q^— ^ ln(iJ^),rf(Q^ y) = Jim Q^-p^\ndDs.an{Q'^ ,y), (52) 



using (^) we find 



Inserting expression ( pB| ) in Eq. ( P^ we obtain the corresponding RG improved expression 
for the propagator amplitude 

D^iQ) , (54) 

here, the factor CRc{y) is determined by the normalization condition (P^). We see that, in 
the weak-coupling limit — (for Nf < 9) expression (0) reproduces the standard RG 
improved solution (|38|) . Furthermore, it has the correct "abelian limit" asdRcijj^y) — > 

«an(x2) ^ ^ 1- 

For ^ > 0.5, the function (Q) satisfies the Kallen-Lehmann analyticity. The corre- 
sponding flavor condition is Nf < 3. For Nf > 3, the unphysical singularities appear 
in the first Riemann sheet. This limitation for Nf turns out to be natural. Indeed, the 



singularities in (|5J) occurred at |^| = 1 where the number of active quarks are just three. 
On the other hand, the flavor condition Nf < 3 seems to be plausible for massless QCD. 
For this reason we cannot reject the solution ([5^ ) using arguments of analyticity |. 



"^The author owes to D.V. Shirkov for drawing his attention to this pecuharity of the solution (|54|). 
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It is interesting to compare the analytic solutions (0), p5| ) and (0), which have 
correct "abelian limit". Only (^) and (54) satisfy an unsubtracted DR. The convenient 
criterion has been formulated in Refs.[Q. This is the principle of minimality for the non- 
perturbative contributions in perturbative (ultraviolet) region. According to this principle 
one can easily verify that the solution (^) is preferable. Indeed, it predicts more rapid 
decrease of the nonperturbative contributions in the ultraviolet region then the solutions 
(^OD , (^^ and (^). However, we cannot accept a final decision for selecting the solu- 
tions. Indeed, the above mention principle is relevant to the full gluon propagator, which 
may include a pure nonperturbative contributions. In the lAA these contributions are 
invisible. 



5 Conclusion 

The RG equation for the QCD two-loop invariant coupling has been solved explicitly. 
We have expressed the solution in terms of the Lambert W function. This allows us to 
understand more clearly the analytical structure of the solution in the complex plane. 
The corresponding analytic coupling has been reconstructed via the dispersion relation. 



It has been demonstrated that the "analyticized" iterative solution (12) is numerically 
close to the "analyticized" exact one ([T7|) . 

We have expressed the invariant (running) coupling of the special RS via the Lambert 
W function. The corresponding analytic invariant coupling is constructed (see (p3D). 
We have shown that the analytic (running) couplings of the two considered schemes are 
numerically close in the IR region. 

The structure of the /5-function has been analyzed in lAA. We have solved the RG 
equation, with nonperturbative model /3-function, giving explicit expression for the in- 
variant coupling as a function of the scale also in terms of the Lambert W function (see 
(|3TD). We have found that the solution is automatically causal. 

The one-loop gluon propagator amplitude of massless QCD is considered in the Landau 
gauge. The RG and analyticity constraints alone are not sufficient to uniquely determine 
the analytic solution for the gluon propagator starting from PT. Therefore, several versions 
of "analyticization" of the gluon propagator are considered. Properties of the obtained 
analytical solutions for the propagator are discussed. 

Finally, we remark that the gluon propagator is central object in the framework based 
on the SDEs |T2| , p!3| , Here, the analytic perturbative solutions can be used to derive 
more complete (nonperturbative) approximants to the gluon propagator. For related ideas 
and applications, see [^ ]. 
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